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Abstract 

Accurate knowledge about the neutron skin thickness IL.Rnp in has far-reaching implica¬ 

tions for different communities of nuclear physics and astrophysics. Yet, the novel Lead Radius 
Experiment (PREX) did not yield stringent constraint on the ARnp recently. We employ a more 
practicable strategy currently to probe the neutron skin thickness of ^^®Pb based on a high linear 
correlation between the ARnp and J—asym, where J and Osym are the symmetry energy (coefficient) 
of nuclear matter at saturation density and of ^^®Pb. An accurate J — Osym thus places a strong 
constraint on the ARnp- Compared with the parity-violating asymmetry Apv in the PREX, the 
reliably experimental information on the J — agym is much more easily available attributed to a 
wealth of measured data on nuclear masses and on decay energies. The density dependence of the 
symmetry energy is also well constrained with the J —Osym- Finally, with a ‘tomoscan’ method, we 
find that one just needs to measure the nucleon densities in ^*^®Pb starting from Rm = 7.61 ± 0.04 
fm to obtain the ARnp in hadron scattering experiments, regardless of its interior profile that is 
hampered by the strong absorption. 
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I. INTRODUCTION 


The nuclear physics overlaps and interacts with astrophysics not only expands its research 
space but also promotes the development of fundamental physics. A great of attention has 
been paid to the equation of state (EOS) of isospin asymmetric nuclear matter in both 
the two helds as the development of radioactive beam facilities and astronomical observa¬ 
tion facilities over the past decade. The symmetry energy that characterizes the isospin 
dependence of the EOS, is a quantity of critical importance due to its many-sided roles in 
nuclear physics [l-7| and astrophysics js-ld]. Although great efforts have been made and 
considerable progresses have been achieved both theoretically and experimentally, its density 
dependence ultimately remains unsolved because of the incomplete knowledge of the nuclear 
force as well as the complexity of many-body systems. Nevertheless, many important and 
leading issues in nuclear astrophysics require the accurate knowledge about it ungently at 
present. 

The symmetry energy S{p) of nuclear matter is usually expanded around saturation 
density po as 

Asym f P Po 


S{p) = J + ^ 
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where J = S'(po) is the symmetry energy at po. The slope parameter L = 3p5S'(p)/5p|po and 
curvature parameter K^yra = 9p^^^*S'/9p^|po characterize the density-dependent behavior of 
the symmetry energy around po. Extensive independent studies have been performed to 


constrain the slope L, but the uncertainty is still large 


15|-ll8|. 


It has been established that the slope parameter L is strongly correlated linearly with the 
neutron skin thickness ARnp of heavy nuclei 


19 


2l[ |. Although the theoretical predictions on 


L and AR^p are extremely diverse, this linear correlation is universal in the realm of widely 
different mean-held models 2^. Accordingly, a measurement of ARnp with a high accu¬ 


racy is of enormous signihcance to constrain the density-dependent behavior of S{p) around 
Pq. Actually, many experimentalists have been concentrating on it with different methods 


including the x-ray cascade of antiprotonic atoms 
ton elastic scattering 


23| , pyg my dipole resonance [2J, [2^ , pro- 


26| . proton inelastic scattering 


27( 1 and electric dipole polarizability 


28|. However, systematic uncertainties associated with various model assumptions are un¬ 


avoidable. The parity-violating electron elastic scattering measurement in the parity radius 
experiment (PREX) at the Jefferson Laboratory combined with the fact that the parity- 


2 
























violating asymmetry A-p\ is strongly correlated with the neutron rms radius, determined 
the A.Rnp to be O-SS^Q^g fm with a large central value compared to other measurements and 
analyses 291. Although it was suggested that ruling out a thick neutron skin in seems 


premature 


3Cll |. in any case, the large uncertainty seems to be not of much help to explore 


the symmetry energy and other interesting issues. In this work, a more practicable strategy 
compared with the PREX at current is introduced to probe the ARnp of ^‘^^Pb together with 


nuclear matter symmetry energy. A new insight into the neutron skin is also provided. 


II. NEUTRON SKIN THICKNESS ARnp PROBED BY THE J - 


The neutron skin thickness of nuclei 


IS given 


as 


AR. 


np 


^5 “ asym(A))A^/^(/ — Jc) — e^Z/(70J)] + Ssw in the nuclear droplet model j^, 32| 

with isospin asymmetry /, nuclear radius R = ^mi a correction R = e^Z/(20JR) 

due to the Coulomb interaction. Z, A are the proton and mass numbers, respectively. Ssw 
is a correction caused by an eventual difference in the surface widths of nucleon density 
prohles. asym(A) is symmetry energy (coefficient) that has been received great interest 
;h the help of it one may obtain some information on the density dependence 


because 

wi 

h t 

of ^(p) 
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35l |. Centelles et al. showed that the neutron skin thickness ARnp correlates 


linearly with J — asym(A) based on different mean-held models, where the symmetry 
energy (coefficient) OgyiRA) is obtained within the asymmetric semi-inhnite nuclear matter 
(ASINM) calculations 3^. In our previous work, instead of using the ASINM calculations, 
the asym(A) was obtained in the framework of the Skyrme energy-density functional 
approach by directly integrating the density functional of the symmetry energy after 


subtracting Coulomb polarization effect without introducing additional assumptions 


33| 


In the present work, the asym(^) of ^°®Pb, marked as agym, is extracted with both the 
Skyrme effective interactions and relativistic effective interaction Lagrangians, and the local 
de nsity approximation is adopted by dropping the negligible non-local terms compared 
As done in Ref. 22j, to prevent eventual biases, we avoid including more than 


to 


two models of the same kind htted by the same authors and protocol and avoid models 
providing a charge radius of ^ospp away from experiment data by more than 1%. 

The calculated neutron skin thickness ARnp of ^‘^^Pb and J — agym with different mean- 


held models are presented in Fig. 1, in which a close dependence of ARnp on J 


^sym 
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FIG. 1: (Color Online) Neutron skin thickness ARnp in 208p^ against the J — asym with different 
nuclear energy-density functionals. 


predicted by the droplet model is displayed. By performing a two-parameter fitting, the 
correlation is given by 


ARnp = (0.0138 ± 0.0003)(J - + (0.0376 ^ 0.0041), 


( 2 ) 


with the correlation coefficient r = 0.989, where AR^p and J — asym are in nnits of fm 


and MeV, respectively. Here the empirical satnration density po = 0.16 fm 


-3 


is nsed 


nniformly. If the symmetry energy is calcnlated at their own saturation densities from 
the mean-field models, the linear correlation vanishes due to the fact that the relativistic 
interactions provide smaller saturation densities compared with the non-relativistic ones. 
The ARnp of ^osp]^ jg found to have a high linear correlation with J — asym as that with the 
slope L (not shown here). It is thus indisputable that the J—agym with a high accuracy places 
a stringent constrain on the ARnp. As the primary advantage, the reliably experimental 
information about the J — agym is much more easily available compared with that about 
the parity-violating asymmetry Apv in the PREX. Recently, the symmetry energy J at 
saturation density po has been well determined to rather narrow regions, in particular, 32.5± 


0.5 MeV from the mass systematics 


371 and 32.10 ± 0.31 MeV from the double differences 


of experimental symmetry energies 38j agreeing with that of the mass systematics. These 
results are very useful in exploring the density-dependent symmetry energy as inputs 
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Here we adopt the union of the two values, i.e. J = 32.4 ± 0.6 MeV, and hence the central 


issue is to determine the symmetry energy Osym of accurately. We extract the mass 

dependent symmetry energy asym(H) = J/(l + ^ l4l| with /3“-decay energies Q/?- 

of heavy odd-H nuclei and with mass differences AH between — 1) and + 1) as our 


42 


43| but with a new input quantity J, and then derive the a, 


sym 


previous calculations 

of The merit of these two approaches is that only the well known Coulomb energy 


survives in Qp- and in AH when determining the unknown Osym, where the and AH 
are all taken from experimental data. Consequently, the agym is extracted to be 22.4 ± 0.4 
MeV accurately, which is quite insensitive to the input quantity J. As a result, the derived 
J — flsym is 10.0 ± 1.0 MeV (solid circle in Fig. 1), which allows us to constrain the neutron 
skin thickness as well as the slope L in our subsequent calculations. 

The neutron skin thickness in jg predicted to be ARnp = 0.176 ± 0.021 fm (solid 

square in Fig. 1), where the estimated error stems from the uncertainties of the J 


■^'sym 


as well as Eq. ([2]). To reach such an error level, the Apv in the PREX should be measured 
at least up to 2% accuracy, which is hardly implemented at present. This fact indicates 


the J — agym is much more effective to probe the ARnp currently. The precise information 
about the ARnp is of fundamental importance and has far-reaching implications in neutron 
star physics, such as the structure, composition and cooling. As an example, a relation of 
pc ~ 0.16 — 0.39AH„p was put forward to describe the relation between the AR ^p o f ^ospp 
and the transition density pc from a solid neutron star crust to the liquid interior where 
the Pc is estimated to be 0.091 ± 0.008 fm“^. The properties of the crust-core transition is 
of crucial importance in understanding of the pulsar glitch 451. 


III. DENSITY DEPENDENCE OF THE SYMMETRY ENERGY PROBED BY 

THE J - agym 

Since the neutron skin thickness ARnp correlates linearly with both the slope L and J — Osym, 
the slope L naturally correlates linearly with the J — agym, which is displayed in Fig. 2(a). 
The linear relation is L = (9.682±0.285)(J —asym) + (—42.694 =f 3.441), where L and J — agym 
are in units of MeV. Imposing the above obtained J — agym, the slope parameter is estimated 
to be L = 54 ± 16 MeV. Recently, the properties of nuclear matter at subsaturation density 
p ^ 0.11fm“^ have attracted considerable attention because it has been shown that the ARnp 
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is uniquely fixed by the slope L{p k. 0.11 fm 4^ and the giant monopole resonance of 


heavy nuclei is constrained by the nuclear matter EOS at this density 


£7|. Fig. 2(b) shows 


that the slope L{p = O.llfm (labeled Lq.h for short) and J 


^sym 


have a higher linear 


dependence Lq.ii = (4.542±0.073)(J—Osym) + (2.140 =f 0.885) with the correlation coefficient 
r = 0.995. Accordingly, the Lq.h is evaluated to be 48 ± 6 MeV. 
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FIG. 2: (Color Online) Correlation of the slope parameter L at densities p = 0.16 fm ^ and 
p = 0.11 fm“^ with the J — Osym- 


The slope L is constrained with the J — Osym in another way for comparison. Centelles 
et al. found that the symmetry energy agym of ^ospp jg approximately equal to the nuclear 
matter symmetry energy S{pa) at a reference density p^ — 0.1 fm“^ 3^. This important 
relation bridges the symmetry energies of nuclear matter and the hnite nucleus. We calculate 
the reference density pA for ^ospp interactions which provide the values of 

J and Osym agreeing with the ones extracted from experimental information, give pa — 0.088 
fm“3 = 0.55po. It should be noted that the Osym does not equal the symmetry energy at 
the mean density of ^ospp as a result of the extremely inhomogeneous isospin asymmetry 
distribution in the nucleus as shown in 33(]. Since the accurate value of the reference density 
Pa is of crucial importance for determining the slope parameter L 42|,l43|], we further examine 
it. Instead of the DDM3Y-shape expression used before 3, 3,Eq. ([1]) is employed directly 
to describe the density dependent symmetry energy to reduce the uncertain factors as far 
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as possible. The ii'sym term that contributes weakly to the symmetry energy nearby po is 


estimated with the relation i^sym = 39 + 5L — 15J 


^1 obtained from the DDM3Y-shape 


expression without loss of accuracy. In terms of J—S{pa) = 10.Oil.0 MeV and pa = 0.55po, 
the slope L at the saturation density po is predicted to be 53 i 10 MeV according to Eq. 
(IT]) , which is in excellent agreement with that from Fig. 2(a). At the density of p = 0.11 
fm“^, the slope To.ii = 49 i 4 MeV, being also particularly consistent with the value of 
48 ± 6 MeV from Fig. 2(b). The consistency of the two approaches not only indicates 
the reliability of the present methods but also further verihes the accuracy of the reference 
density pa = 0.55po. As an important conclusion, the Osym = S{p = 0.55po) — 22.4 MeV 
will be a very useful reference to calibrate the effective interactions in nuclear energy density 
functionals. 

With the obtained Lo.n and L values, the curvature parameter is evaluated to be 
A'sym = —152 ± 70 MeV. Currently, the symmetry energy at suprasaturation densities is 
extremely controversial. It was indicated that the three bulk parameters J, L and iCgym 
well characterize the symm etry energy at densities up to ~ 2po while higher order terms 
contribute negligibly small [49[. If true, the symmetry energy S'(p) at high densities up to 
~ 2po turns out to be not stiff, as shown in Fig. 3. The symmetry energy at 2po is estimated 
to be S'(2po) = 42±10 MeV. In short, to characterize the symmetry energy at high densities, 
the accurate knowledge about its density dependence at the saturation density is crucial. 



FIG. 3: (Color Online) Density dependent symmetry energy at high densities. 
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IV. FURTHER EXPLORATION ON THE MEASUREMENT OF THE AR, 


np 


Based on the above discussions on the neutron skin thickness ARnp and symmetry energy, we 
make an exploration on the measurement of the ARnp in To grasp richer information 

on the ARnp, we formulate it as an integral of a distribution function 

^Rnp = V<^1> - f{r)dr, (3) 

where /(r) = ~ ^) / ^a/< > + ^/< rp is dehned as the radial distribution 

function which is actually determined by the nucleon densities and reflects the detailed 
information about the neutron skin. \/< t‘^> + a/< ^p > — H-l f^i changes by less than 
3% in the mean held model calculations, and can be taken as a known value. Fig. 4(a) 
illustrates the distribution function /(r) in ^°®Pb as a function of distance r generated by 
the SLy5 interaction as an example. It is a misleading idea to consider the neutron skin 
merely originating from the nuclear surface. The area enclosed by the x-axis and the curve 
/(r) (colored regions) is exactly the neutron skin thickness Ai?„p. We name this new method 
that dissects the Ai?„p with a distribution function as ‘tomoscan’ picturesquely here. As a 
new concept in nuclear physics, it could also be used to analyze other intriguing issues, such 
as the halo structure in exotic nuclei. The region of r < i?o contributes negatively while 
that of r > i?o contributes positively to the ARnp. Thus, there exists a distance Rm below 
which (0 < r < R^) the contributions (red shaded regions) cancel each other out, and hence 
the ARnp can be calculated by the neutron and proton density distributions just starting 
from Rm (blue hlled region). 

The calculated values of Rm with different interactions are marked in Fig. 4(b). The 
Rm is found to be model dependent, which should be further constrained. The interactions 
generating smaller (larger) ARnp tend to yield slightly larger (smaller) Rm- As we mentioned 
above, one important conclusion of this work is that the Osym = = 0.55po) — 22.4 

MeV (along with J ~ 32.4 MeV) serve important calibrations for effective interactions in 
nuclear energy density functionals. Thus we use those constraint conditions to hlter those 
interactions. The eligible interactions give Rm = 7.61 ± 0.04 fm (colored solid symbols), 
where the error bar of ±0.04 fm just leads to an uncertainty of the ARnp by about ±0.005 
fm. The error bar of ±0.005 fm for the ARnp is so small that the obtained Rm value 
should not be regarded as model dependent any more. This result leads to an intriguing 
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FIG. 4: (Color Online) (a) Radial distribution function /(r) of the neutron skin thickness in 
^^®Pb. The contributions from the two parts in the red shaded regions cancel each other out. The 
area under the curve of /(r) starting from (blue filled region) is equal to the neutron skin 
thickness ARnp- (b) Calculated Rm values with different energy density functionals. The colored 
solid symbols are from the interactions generating the reference density pA — 0.55po, Osym — 22.4 
MeV and J ~ 32.4 MeV. (c) Error accumulation of the ARnp measurement in hadron scattering 
experiments as a function of distance r, where the nucleonic density distributions are from Tables 


III and IV in Ref. 


m 


conclusion: one just needs to measure the rather dilute matter located in the nuclear surface 
to determine the neutron skin thickness of namely, only measures the nucleon densities 

from r = = 7.61 ± 0.04 fm to about r = 12 fm. Thus, the measurement of the 

ARnp would be substantially simplihed in hadron scattering experiments which have been 
hampered by the strong absorption in the nuclear interior. We stress that, contrary to the 
usual understanding, the nuclear surface properties are in fact not well constrained by the 
nuclear mean-field models obtained by fitting nuclear masses and charge radii. For instance, 
both the SLy5 and NL3 interactions give Rm = 7.62 fm, but they provide a substantial 
difference in the ARnp amounting to 0.12 fm. In other words, it is exactly the ambiguity 
of the nuclear surface profile that leads to the large uncertainty of the ARnp, because the 
radial distribution function /(r) relies on the fourth power of distance r according to Eq. 
(3), causing a drastic amplification of the error as r increases. Fig. 4(c) illustrates the error 
accumulation of the ARnp in hadron scattering experiments for different regions, which is 
obtained by analyzing the data in Ref. 2^ combined with the ‘tomoscan’ method. The 
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error accumulation at distance r is defined as the error generated by the region from the 
nuclear center to r. It indicates that the error also primarily originates from the surface 
structure. Therefore, the surface prohles must receive particular attention and be measured 
with a much higher accuracy. 


V. SUMMARY 

We have developed alternative methods in the present study to explore the neutron skin 
thickness of ^°®Pb and density dependence of symmetry energy. The main conclusions 

are summarized as follows, i) We have established a high linear correlation between the 
and J —Ogym on the basis of widely different nuclear energy-density functionals. Accordingly, 
an accurate J — Osym value sets a signihcant constrain on the ARnp, which turns out to 
be a much more effective probe than the parity-violating asymmetry Apv in the current 
PREX. ii) The symmetry energy (coefficient) Osym of ^°®Pb was extracted accurately with 
the experimental /9“-decay energies of heavy odd-A nuclei and with the experimental mass 
differences. Given that the symmetry energy J has been well determined recently, the ARnp 
in 208pb 

was thus predicted to be 0.176 ± 0.021 fm robustly. This conclusion would be 
signihcantly meaningful to discriminate between the models and predictions relevant for the 
description of nuclear properties and neutron stars, iii) With the above derived J — Ogym, 
the values of the slope L of the symmetry energy at the densities of p = 0.16 fm“^ and 
p = 0.11 fm“^ which are of great concern, are predicted to be 54 ± 16 MeV and 48 ± 6 MeV 
respectively. These results, together with the ARnp of ^°®Pb, can be applied to explore some 
intriguing problems in nuclear astrophysics. In particular, the derived Osym and S{pa) serve 
as important calibrations for a reliable construction of new effective interactions in nuclear 
many-body models, iv) The symmetry energy at suprasaturation densities up to ~ 2po was 
predicted to be not stiff, v) With the hrstly proposed ‘tomoscan’ method, we concluded 
that to obtain the ARnp one needs to only measure the nucleon densities in ^°®Pb from 
Rm = 7.61 ± 0.04 fm as the densities in the range of r < Rm have no contribution to the 
ARnp. Thus, the measurement on the ARnp is signihcantly simplihed in hadron scattering 
experiments which have been hampered by the strong absorption in the nuclear interior. 
Incidentally, the ‘tomoscan’ method could be employed to analyze the halo structure in exotic 
nuclei, vi) It has been widely believed that the nuclear surface structure is well constrained 
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in nuclear energy-density functionals and in experimental measurements. However, within 
the ‘tomoscan’ concept, we have showed that it is not true but a complete illusion. To grasp 
the ARnp, one must especially concentrate on the dilute matter located in nuclear surface 
which results in the dominant uncertainty. 
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